Green Function theory vs. Quantum Monte Carlo Calculation for thin magnetic films 
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In this work we compare numerically exact Quantum Monte Carlo (QMC) calculations and Green 
function theory (GFT) calculations of thin ferromagnetic films including second order anisotropies. 
Thereby we concentrate on easy plane systems, i.e. systems for which the anisotropy favors a 
magnetization parallel to the film plane. We discuss these systems in perpendicular external field, i.e. 
B parallel to the film normal. GFT results are in good agreement with QMC for high enough fields 
and temperatures. Below a critical field or a critical temperature no collinear stable magnetization 
exists in GFT. On the other hand QMC gives finite magnetization even below those critical values. 
This indicates that there occurs a transition from non-collinear to collinear configurations with 
increasing field or temperature. For slightly tilted external fields a rotation of magnetization from 
out-of-plane to in-plane orientation is found with decreasing temperature. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.70.Ak, 75.30.Gw 
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I. INTRODUCTION 

The fast development of technological applications 
based on magnetic systems in the last years, e.g. 
magnetic data storage devices, causes a high interest 
in thin magnetic films. One precondition for the tech- 
nological development is the investigation of magnetic 
anisotropies and spin reorientation transitions connected 
therewith. Those reorientation transitions can occur 
from out-of-plane to in-plane or vice versa for increasing 
film thickness c&, temperature T2£&&&2& 1 or external 
field B . 

Quantum Monte Carlo (QMC) calculations give the 
possibility to compare numerically exact results with 
analytical approximations. In Ref. HI the authors inves- 
tigated a ferromagnetic monolayer including positive 
second order anisotropy (easy axis perpendicular to the 
film plane). They discuss the temperature dependence 
of the magnetization (S Z )(T) as well as field induced 
reorientation transitions from out-of-plane to in-plane 
and compare the QMC results with Green function 
theory (GFT). They found good agreement in the case 
of applied external field in the easy direction (here 
z-axis). However, their GFT fails for external field 
applied in arbitrary direction, especially in the hard 
direction (within the film plane). As shown in Ref. 
for getting closer to the QMC results for magnetic field 
induced reorientation from out-of-plane to in-plane a 
more careful treatment of the local anisotropy terms is 
needed. In Refs. [i3Tl1ll^ll3l a decoupling scheme was 
presented which yields excellent agreement with QMC 
results for out-of-plane systems. 

The availability of theories such as GFT and their 
check against state-of-the-art numerical algorithms 
is highly desirable because of the size limitations of 
systems where QMC can be performed. On the other 
hand the extension of GFT from a monolayer (where it 



can be compared to QMC as in the the present work) 
to multilayer systems is a straightforward task without 
further approximations 1 ^. 

Up to now, to our knowledge, there is no comparison 
between QMC and approximative theories for easy-plane 
systems and it is not obvious that the theory presented 
in Refs. [IMIIIGJ can reproduce the QMC results for 
in-plane systems as accurately as for the out-of-plane 
case. In contrast to the easy-axis case where a certain 
direction is preferred by the single ion anisotropy in 
easy-plane systems the full xy-plane is favored and no 
particular direction is distinguished within the plane. A 
magnetic field applied perpendicular to the plane does 
not destroy the xy-symmetry. 

For systems exhibiting this kind of symmetry it was 
shown in a classical treatment that for external fields 
smaller than a critical field < B < B cr n (B || z) 
stable vortices, i.e. a non-collinear arrangement of 
spins, can exis t 15 ! 16 ' 17 : 18 ' 19 . These vortices can undergo 
a Berezinskii-Kosterlitz-Thouless (BKT) transition 14 . 
Depending on the strength of the anisotropy there 
might be vortices with or without a finite z-component 
of magnetization 15 . In the small anisotropy case (which 
is considered in this work, \K%\ < 0.1J) there is a 
finite out-of-plane component and for zero field the 
two possible directions of magnetization (±z) are en- 
ergetically degenerate. For increasing magnetic field in 
z-direction the vortices antiparallel to the field become 
more and more unstable (heavy vortices). However the 
so called light vortices (parallel to the field) are stable 
up to a critical field B z — B cr u and contribute a finite 
z-component to the net-magnetization of the considered 
system^ 9 .. 

The vortices in connection with a finite z-component 
of the net-magnetization emerge because of two reasons: 
first the competition between the anisotropy (favoring 
a orientation of the magnetization within the sy-planc) 
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and the external field (favoring a perpendicular magne- 
tization), and second: the xy-symmetry of the system, 
which does not allow for a rotated homogeneous phase. 

In this paper we investigate both aspects, i.e. the 
field vs. anisotropy competition as well as the symmetry 
properties in detail for a quantum mechanical system. 
We will compare the results of QMC and GFT calcula- 
tions. 

As explained in more detail below, the QMC al- 
gorithm used here allows only for an external field 
applied in z-direction. Thus the xy-symmetry can not 
be broken and no comparison between xy-symmetric 
and asymmetric systems is possible. We will use GFT 
to clarify the influence of this symmetry breaking on 
the homogeneous phase. On the other hand, the GFT 
used here is by ansatz limited to the homogeneous 
phase. Therefore it can not describe a non-collinear 
(e.g, vortex-) magnetic phase, which is expected for 
B || z and small field strengths. The breakdown of 
magnetization in GFT as well as an exposed maximum 
in the magnetization in QMC at certain critical values 
of the external field or temperature gives, however, a 
clear fingerprint of non-collinear configurations, at least 
if there is no meta-stable homogeneous phase. Below 
these critical values there will be a finite z-component 
in QMC and a vanishing magnetization in GFT. 
For parameters, where both theories are applicable, 
QMC serves as a test for the approximations needed in 
GFT. 

In this work we find indications for non-collinear 
spin configurations below a critical field or temperature 
for B || z by comparing results of QMC and GFT as 
explained in the last clause. Above the critical field 
we obtain good agreement between QMC and GFT 
results. Breaking the xy-symmetry by adding a small 
x-component to the external field yields a stable collinear 
solution in GFT. The z-component of the magnetization 
in this case is in good agreement with the QMC results 
calculated with untilted field. Thus we can conclude 
that except for the restriction to collinear magnetic 
states GFT describes the competition between external 
field and anisotropy quite well. 

The paper is organized as follows: First we ex- 
plain the basics of the GFT and the QMC calculations. 
Then we apply both approaches to easy-plane systems 
in external magnetic fields and report the results of our 
calculations. 



II. THEORY 

A. Green Function Theory 

In the following we present our theoretical approach us- 
ing Green function theory. The focus of this work lies on 
the translational invariant system of a two-dimensional 



monolayer. Therefore the following Hamiltonian is used: 
H = -^JyS^ -B^S,-^^^) 2 . (1) 

ij i i 

The first term describes the Heisenberg coupling Jjj be- 
tween spins Si and Sj located at sites i and j. The second 
term contains an external magnetic field B in arbitrary 
direction (the Lande factor gj and the Bohr magneton 
Hb are absorbed in B). The third term represents second 
order lattice anisotropy due to spin-orbit coupling. S z ,i 
is the z-component of S; (the z-axis of the coordinate 
system is oriented perpendicular to the film-plane). The 
lattice anisotropy favours in-plane (K2 < 0) or out-of- 
plane (K2 > 0) orientation. Our Hamiltonian is similar 
to that used in Refs. [T^IUIIIIilH! for the investigation 
of the magnetic anisotropy and the field induced reori- 
entation transition. To simplify calculations we consider 
nearest neighbor coupling only 

J = { J W'O') n - n - (2) 
v \ otherwise. ^ ' 

The main idea of the special treatment presented in 
Refs . [IollTTIIT^fl3l is that, before any decoupling is applied, 
the coordinate system £ is rotated to a new system E' 
where the new z'-axis is parallel to the magnetization im- 
plying a collinear alignment of all spins within the layer. 
Then a combination of Random Phase approximation 
(RPA}2i for the nonlocal terms in Eq. |T]) (Heisenberg 
exchange interaction term) and Anderson-Callen approx- 
imation (AC)2£ for the local lattice anisotropy term is 
applied in the rotated system. After application of the 
approximation one gets an effective anisotropy 

K eff {T) = 2K 2 (l-^2 ( S ( S + !) - (#») < S ») ( 3 ) 

where (S zl ) is the norm of the magnetization and S is the 
spin quantum number, that we have chosen to be S = 1 
in all our calculations. 

As shown in comparison with an exact treatment of the 
local anisotropy term in Ref. [26| this approximation still 
holds up to anisotropy strengths K 2 ~ 1/2 J. Therefore 
we restrict ourselves in the following to small anisotropics 
(K2 < 0.1 J) as found in most real materials^ 3 .. For a 
magnetic field applied in the xz-planc (B = (B x ,0, B z )) 
our theory gives a condition for the polar angle 9 of the 
magnetization: 

sin 9B Z - cos 6B X + K ef / sin 9 cos 9 = (4) 

The uniform magnon energies (q = 0) which dominate 
the physical behavior of the magnetic system can easily 
be extracted from the theor y 12 ' 13 : 

El =0 = (cos 9B Z + sin6B x + K eff (cos 2 9 - sin 2 9)) ■ 
(cos 9B Z + sin 9B X + K ef f cos 2 9) (5) 
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This result coincides with the spin-wave result^ if one 
replaces (S zl ) by the spin quantum number S and K e ff 
by the bare anisotropy constant K<z in Eq. (|5|). For an 
easy-plane system (K e ff < 0) with external field B in 
^-direction the polar angle of the magnetization^ is 
given by: 



cos 9 



-B/K eff (T) 
1 



for B < \K eff (T)\ 
otherwise 



(6) 



By inserting Eq.© into Eq.© one immediately gets: 



B 





Keff(T) 



B<\K eff (T)\ 
otherwise. 



(J) 



For gapless magnon energies -E q= n = the magnon oc- 
cupation number <fi diverges (cf> — > oo) in film systems 
with ferromagnetic coupling J > and the magnetiza- 
tion becomes zero {S z/ ) — in the collinear phase. This 
can be seen by following an argument of Bloch-2, already 
given in 1930. Since the spin wave dispersion is E w q 2 
in the vicinity of q = the spin-wave density of states 
N(E) is independent of E for a two-dimensional system 
for E close to zero. The excitation of spin-waves at finite 
temperature leads to a variation of the magnetization of 
the order: 



Am(T) 



N(E)dE 



, exp(E/k B T) - 1) 

k B T / 

Jo 



exp(a;) — 1 



(8) 



Since the integral in Eq. ([5]) diverges for T ^ and exited 
spin-waves lead to a reduction of the magnetization one 
can conclude that the magnetization should be zero at 
finite temperature. However for an infinitesimally small 
contribution of the external field parallel to the plane, 
i.e. B x ^ 0, a finite gap in the excitation spectrum 
at q = opens. This can be seen in FigQ] where the 
uniform magnon modes -E q= o {B) are shown for different 
orientations 6b, where 9 b is the polar angle of the ex- 
ternal field. The integral f8j is now finite and a stable 
finite magnetization in the collinear phase having a well 
defined orientation in the xz-plane is possible. 

Let us now come back to the case where the applied 
field is aligned in z-direction. It can be seen from Eq. 
(J7|) that for external field B (B \ \ z) larger than a critical 
field B > B crit given by: 



B cnt = \K eff (T,B) 



(9) 



a stable collinear solution exists. Since K e ff(T) is a de- 
creasing function of temperature T a transition from non- 
collinear to collinear phase with increasing temperature 
is possible. In Fig. f2] we show the normalized critical 
field © B crit / K 2 as a function of temperature T. For 
a constant magnetic field B (B \\ z) at a temperature 
Ti with B < K e f f (Ti , B) no stable collinear phase ex- 
ist. Then by increasing the temperature up to T 2 the ef- 
fective anisotropy K e ff is sufficiently reduced such that 
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FIG. 1: The energies of the uniform magnon mode Eq = o(B) 
for different polar angles 9b of the external field. -B q= o is 
zero below B/J ~ 0.03 for 9b =0°. The prefactors gj\iB 
and fcs are absorbed in B and T respectively. The latter are 
given in units of the nearest neighbor Heisenberg coupling J. 
Parameters: 5 = 1, K/J = -0.03 and T/J = 10" 4 . 
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FIG. 2: The normalized critical field Bcrit /Ki as a function 
of temperature. Parameters: 5=1. 



B > K e ff(T2,B), and the collinear phase becomes sta- 
ble. Before we come to the results let us briefly sketch 
the main aspects of the QMC. 



B. QMC 

In the last section we gave a short description of the 
theory used to treat a system described by a Hamilto- 
nian of form |T]). This theory applies to the thermody- 
namic limit (films of infinte size) but contains certain 
approximations. Additionally the GFT is restricted to 
ordered phases with a collinear alignment of all spins. 
Therefore it would be very useful to have exact results at 
hand to crosscheck the predictions of GFT. A Quantum 
Monte Carlo method, particularly well suited for spin 
systems, is the stochastic series expansion (SSE) with 
directed loop update. We will sketch this method here 
only briefly as detailed descriptions can be already found 
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elsewhe re 28 ' 29 ^. 

Our starting point is the series expansion of the parti- 
tion function 

Z = Tre-f" = ]T E ^-<«l(-^ n l«> (10) 

n— a 

where H denotes the Hamiltonian, {\a)} are basis vectors 
of a proper Hilbert space and /3 is the inverse tempera- 
ture. The Hamiltonian is then rewritten in terms of bond 
Hamiltonians: 

M 

H=-jJ2H b (11) 

b=l 

where Hb can be further decomposed into a diagonal and 
an off-diagonal part: 

H D ,b = C + SfaSjQ) +bb[Sfy + Sj^] (12) 
+k 2b [(S* (b) ) 2 + (S* (b) f] 

H ,b = 2^ s i(b) s m + S i(b) s t(b)} ( 13 ) 

Here we have renormalized the anisotropy constant kib 
and the magnetic field bb in such a way that pip coincides 
with ([1]). i(b) and denotes the lattice sites connected 
by the bond b and the additional constant C in H^.b will 
be chosen such that all matrix elements of this term be- 
come positive, a condition necessary to interpret them as 
probabilities. Note that for a finite system at finite tem- 
perature the power series of the partition function can 
be truncated at a finite cutoff length A without intro- 
ducing any systematic error in practical computations 29 . 
Therefore reinserting (|11[) into (|10[) and rewriting the re- 
sult yields: 

A v^v^ tf n (A-n)! ; ,„ , v 
Z = EEE A! ( a \ S cM- ( 14 ) 

n=0 Sc A a 

Here Sc A denotes a product of operators (operator 
string) consisting of n non-unity operators and (A — n) 
unity operators Hq = Id which were inserted to get op- 
erator strings of equal length A. 

In fact it is impossible to evaluate all operator strings in 
(fT4")l . The SSE-QMC replaces such an evaluation there- 
fore by importance sampling over the strings according to 
their relative weight. Hence an efficient scheme for gen- 
crating new operator strings is needed. In the directed 
loop version of the SSE this is done by dividing the up- 
date into two parts. In a first step a diagonal update is 
performed by traversing the operator string and replac- 
ing some unity operators by diagonal bond operators and 
vice versa (the probabilities for both substitutions have 
to fulfill the detailed balance criterion). Then the loop 
update follows in which new non-diagonal bond opera- 
tors can appear in the operator string. For details of the 
update procedure we refer the interested reader to the 




T/J 

FIG. 3: Magnetization vs. temperature for an out-of-plane 
easy-axis system (K 2 > 0). Straight line: GFT (RPA+AC) 
result; symbols: QMC results for different system sizes N 2 . 
Parameters: S = 1, B/J = 0.01 (B \\ z) and K 2 /J = 0.01. 

according literatu re 28 ! 29 ' 30 . 

A full implementation of the SSE with directed loop 
update which we have used for all QMC calculations 
in this work can be found in the ALPS projec t 30 ' 31 . 
Since the SSE-QMC used by us is implemented in z- 
representation (spin quantization axis along z-axis) in- 
plane correlation functions e.g. the in-plane magnetiza- 
tion are not accessible. Further B \ \ z is the only possible 
field direction in the used QMC implementation because 
a traverse field (in-plane field component) would lead to 
non-closing loops (see Ref. 0) . 



III. RESULTS 

As mentioned in Sec. Ill Al the results for the in-plane 
systems are very sensitive to the effective anisotropy 
K e ff(T). This sensitivity of the anisotropy is less pro- 
nounced for out-of-plane systems (K 2 > 0) since the ap- 
plied field B (B\\z) and the intrinsic easy axis are par- 
allel. In order to test our decoupling scheme (RPA+AC) 
we first compare GFT and QMC for an out-of-plane 
system^ 

In Fig. [3] the magnetization (S z ) as a function of tem- 
perature T is shown. The straight line belongs to the 
GFT whereas the symbols show the result of the QMC 
for different system sizes. Let us first comment on finite 
size effects in the QMC results. 

It can be seen in Fig. [3] that the QMC results converge 
for increasing system size A^ 2 (for N x N square lattice) . 
Indeed for A^ > 32 the QMC results are unbiased by finite 
size effects and resulting magnetization curves are almost 
equal for increasing A^ > 32. Note that we have omitted 
error bars in the figures showing QMC results because 
the relative errors are of the order 1CP 4 . 

We now compare the GFT with the QMC results 
(JV = 64). For low temperatures (T/J < 0.5) we ob- 
tain excellent quantitative agreement. This is plausible 
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FIG. 4: z-component of magnetization as a function of ex- 
ternal magnetic field for fixed temperature T/J — 0.4. In 
contrary to the GFT the magnetization obtained by QMC 
remains finite for all fields. The QMC results are converged 
for N > 64. Parameters: S = 1, K%/J = -0.06 and 6 B = 0°. 



FIG. 6: z-component of magnetization vs. external field for 
different temperatures T/J and fixed system size N 2 (N = 
128). Solid lines: GFT (Ob = 0.5°), dashed lines: GFT (Ob = 
0°) other parameters as in Fig. [3] 




FIG. 5: z-component of magnetization vs. external field for 
T/J = 0.4 with slightly tilted field (0 B = 0.5°)in the GFT 
result (solid line). The dotted line shows GFT result for (Ob = 
0°). Other parameters as in Fig. [4] 

because in this region the GFT result coincides with the 
result of the spin-wave theory which is known to be re- 
liable (exact for T = 0) for low temperatures. For the 
intermediate region T/J = 0.5. .1 the RPA slightly un- 
derestimates the magnetization which was also found in 
Ref. [9|. The opposite is the case in the region near the 
extrapolated Curie temperature Tc^, where the magne- 
tization is overestimated. The reason is the presence of 
longitudinal fluctuations, which play an important role 
in this region and it is well known that the RPA fails to 
treat them properly. 

We consider now the case of in-plane systems (Ki < 0) 
and applied field in the hard direction (B \ \ z). As al- 
ready mentioned there is no 'collinear' magnetization in 
the GFT for B z < \K eff (T)\. In Fig. g] the z-component 
of the magnetization is shown as a function of the ex- 
ternal field B for a constant temperature T/J = 0.4. 



As in Fig. [3] we see that the QMC results for N > 64 
are almost converged and the finite size of the calculated 
system in QMC should not influence the results anymore. 
The dotted line marks a critical field B crit . For magnetic 
fields larger than the critical one B > B cr u we obtain 
good agreement between QMC and GFT results. Below 
the critical field B < B cr n GFT does not yield a stable 
homogeneous magnetization. However the QMC results 
show that there is a finite z-component of the magneti- 
zation in the considered system for < B < B cr i t . 

In order to compare QMC with GFT results we have 
tilted the magnetic field by 0b = 0.5° which corresponds 
to B x < 10~ 2 B Z in the GFT. As explained before any 
symmetry breaking field B x ^ leads to a stable homo- 
geneous magnetization with well-defined orientation in 
the a;z-plane. However such a small contribution of the 
external field within the plane should hardly influence the 
z-component of the magnetization. This is confirmed by 
Fig. |5] where we show QMC results (N — 128, 0b = 0°) 
as well as the corresponding GFT results with Ob = 0° 
and 0b — 0.5°. As expected for \B\ > B cr u the two so- 
lutions in the GFT are nearly the same and agree well 
with QMC. Below the critical field only the solution with 
the slightly tilted field yields a stable homogeneous mag- 
netization and its z-component compares well with the 
QMC result in the untilted case. 

The above results can be interpreted within a semi- 
classical picture of non-collinear vortex configurations 
which are stable below a critical field B cr n in z-direction 
and contribute a finite z-component to the magnetiza- 
tion in case of an applied field^ Despite the lack of di- 
rect, quantitative access to such states (or correspond- 
ing physical in-plane observables) within the QMC al- 
gorithm they are included in principle and one can ob- 
serve their consequences, namely a finite z-component of 
the magnetization below the critical GFT field. On the 
other hand GFT can only describe homogeneous collinear 
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FIG. 7: The z-component of magnetization as function of 
temperature for a fixed external field. Below a critical tem- 
perature T cr it there is a breakdown of magnetization in GFT 
where is no in QMC. Parameters: B/J = 0.03, 5 = 1, 
K 2 /J = —0.06. 



configurations of spins therefore showing a breakdown of 
magnetization. However by applying a small field in in- 
direction the xy-symmetry is broken and the spins ro- 
tate in the field direction (the vortices vanish) and the 
collinear phase is retrieved. Our results corroborate this 
interpretation based on the classical picture. Let us em- 
phasize that both, GFT for slightly tilted field and QMC 
for B || z, describe the competition between the external 
field (which favors magnetization parallel to z) and the 
anisotropy favoring in-plane magnetization. Comparing 
the z-components of the magnetization for both cases, 
one can conclude that the ratio of the competing forces 
are comparable for QMC and GFT. This indicates that 
this competition is correctly taken into account in GFT. 

In Fig. [B]the same field dependence of the z-component 
of the magnetization is shown for different temperatures. 
We have plotted the result for the tilted field in case of 
GFT, the point of breakdown in the untilted case is in- 
dicated by the dotted line. It can be seen that for higher 
temperatures no breakdown of collinear magnetization 
occurs, meaning that the condition for the critical field 
(B < \K eff (T, B)\) is never fulfilled in this case. The dis- 
crepancies at intermediate temperatures (T — 0.9. .1.2) 
are due to the RPA decoupling in the GFT as was dis- 
cussed already. 

In Figs. [3 [8] and [9] the z-component of the magnetiza- 
tion is plotted as a function of temperature obtained by 
GFT (straight line RPA+AC) as well as QMC (symbols) 
for different system sizes and a constant applied magnetic 
field. 

Let us first discuss the qualitative behavior of the mag- 
netization as a function of temperature which is found in 
all three figures. For high T (T ^> T cr u) the magneti- 
zation is reduced by thermal fluctuations (where the tail 
of the curve above T/J « 1.5 is due to the applied ex- 
ternal field). In the vicinity of T cr n, T — T cr n — > + , a 
competition between two effects sets in and has a pro- 



FIG. 8: Same situation as in Fig. [7]for K 2 /J = -0.04 (other 
parameters as in Fig[7]). The result for a two layer film treated 
by GFT is plotted also (dashed-dotted line). 



nounced influence on the magnetization. On the one 
side the effective anisotropy acts against the external 
field (B eff = B - \K eff (T)\, (B \\ z)). The effective 
anisotropy K e ff(T) is reduced with increasing temper- 
ature T and thus the effective field B e ff increases with 
T. This effect tends to enhance the magnetization with 
T. On the other side thermal fluctuations suppress the 
magnetization with increasing T. The flattening of the 
magnetization curve near T cr a is a result of this com- 
petition. For low temperatures T < T cr n the effective 
anisotropy in the GFT cannot be overcome by the exter- 
nal field (B < \K eff (T)\, (B\\z)). Therefore the collinear 
magnetization in our approximation vanishes due to the 
mentioned gapless excitations, in contrast to QMC which 
yields again a finite magnetization because non-collincar 
states are taken into account as discussed above. The 
reduction of the z-component of magnetization in QMC 
below T cr it can be pictured classically as the spins being 
in a non-collinear phase with an angle 9 with respect to 
the z-axis. Since in general anisotropy effects (which fa- 
vor in-plane magnetization) increase when temperature is 
lowered the z-component of the magnetization decreases. 

Now we discuss the three figures in detail. In Fig. [7] 
we have plotted QMC results for different system size 
showing again that these are well converged for N > 64. 
Thus we conclude that the striking difference between 
GFT and QMC is not a mere finite size effect. The 
breakdown of magnetization in GFT occurs at a critical 
temperature T cr i t /J = 0.5 whereas no such breakdown 
exists in QMC. However the exposed maximum of the 
magnetization in QMC lies near the breakdown point. 
The differences between QMC and GFT in the tempera- 
ture range T/J ps 0.3 . . . 1.3 are due to the decoupling of 
the exchange and anisotropy term in GFT as also seen 
in Fig. [3] It is worth mentioning that the value of the 
z-component of the magnetization is nearly the same at 
the breakdown point in GFT and the maximum in the 
QMC. Thus we have the result that although GFT can- 
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FIG. 9: Same situation as in Fig. [7]for K 2 /J = -0.01, B/J = 
0.005 and slightly tilted field (6 B = 0.5°) for the GFT results. 



not describe the non-collinear phase by ansatz its break- 
down coincides rather well with the onset of this phase, 
which we attribute to the maximum of the QMC curve. 
FigJS] shows the same situation for a different anisotropy 
constant K 2 — —0.04. The critical temperature is lower 
than in Fig[7] since the ratio B z /K 2 becomes larger. The 
tilted field case is also shown for the GFT results. Again 
the qualitative agreement of the ^-component of magne- 
tization with QMC is good. To confirm this point we 
have plotted the temperature dependence for an other 
set of parameters in Fig. [9] There is as good qualitative 
agreement of the two approaches. Additionally one gets 
a finite component in x-direction in GFT which is also 
plotted in the figure. The two effects of the external field 
vs. anisotropy competition are nicely to be seen: a non- 
collinear state for B | | z (z-component only in QMC but 
not in GFT) and rotation of magnetization for slightly 
tilted external field (seen only in GFT). The ratio of the 
competing forces agree well again in both treatments. 

In Fig. [7] we have plotted the results of a different 
decoupling scheme of the anisotropy terms (namely a 
mean field decoupling, dashed line in Fig. [7]) . Although 
the overall characteristic resembles the RPA+AC result 
(breakdown of magnetization) the mean field results dif- 
fer extremely from the QMC for a large range of tempera- 
ture and underestimates the magnetization. This demon- 
strates the reliability of the Anderson-Callen treatment of 
the local anisotropy terms presented in Refs. Il^llllfl^fl3l . 

The extension of the GFT method to multi-layer films 
is straightforward*^ We have also included results for a 
two-layer film in Fig. [5] for the same parameters as in 
the monolayer case. One finds that for a double layer 
magnetism is stabilized, which can be attributed to the 
increased coordination number and thus higher exchange 
energy. Just like for a monolayer, one observes a break- 
down of collinear magnetization at some critical temper- 
ature. This is due to the fact that the same reasoning 
regarding the vanishing excitation gap also applies for 
multilayer (slab) systems 3 - . The effective anisotropy per 



layer is essentially the same as for a single layer, thus the 
critical (S^) -value (magnetization at critical field B cr a) is 
practically the same. The critical temperature is higher 
than that of a monolayer due to the increased magnetic 
stiffness of the double layer. 

IV. SUMMARY AND CONCLUSIONS 

Using GFT and QMC calculations we studied easy- 
plane systems as well as easy-axis systems with an exter- 
nal field applied perpendicularly to the film. The GFT 
treatment of the Hamiltonian Eq. (Q]) consists of a RPA- 
decoupling for the nonlocal terms and an AC-decoupling 
for the local terms performed in a rotated frame, where 
the new z'-axis is parallel to the magnetization. For the 
QMC calculations we have used the stochastic series ex- 
pansion (SSE) with directed loop updates, which is well 
suited for spin-systems. 

We have calculated the magnetization as a function of 
the external field as well as temperature. We found a 
critical field and critical temperature respectively below 
which is no magnetization in GFT whereas there is one in 
QMC. By tilting the field slightly in GFT so that it has 
a small component in x-direction we get a stable magne- 
tization even below the critical field or temperature. The 
z-component of the magnetization in this case coincides 
well with the z-component obtained by QMC for the un- 
tilted field confirming that GFT and QMC agree well in 
the description of the external field vs. anisotropy com- 
petition. However, this comparison can be only some- 
what indirect, since QMC has access to the non-collinear 
(B 1 1 z) state only, while GFT is limited to collinear fer- 
romagnetic states (rotated homogeneous magnetization) 
found for slightly tilted external fields. 

For parameters that are accessible by both QMC and 
GFT (B\\z;B> B cHt (T)) QMC and GFT are in good 
agreement. Thus one can conclude that the GFT is ap- 
plicable to the homogeneous phases of systems described 
by Eq. ([I]) and can be used also for system configurations 
not accessible by QMC due to too large system size as 
e.g. multilayer systems. 

It would be an interesting task for a succeeding work 
to extend the GFT in order to get a deeper insight into 
the non-collinear configurations also. 
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APPENDIX A: MAGNETIZATION ANGLE 

Here we will discuss the second mathematical solution 
which occurs besides Eq. [5J For an external field in 
the z-direction the angle dependent part of the free en- 
ergy including second order anisotropy can be expanded 
asi^: 



Therefore one gets two solutions for in-plane systems 
(K 2 < 0). For sin 8 ^ one gets immediately the so- 
lution of Eq. [6] if 2K 2 /M Z = K e ff holds. This is the 
stable solution. The trivial (second) solution sin 8 = is 
unstable for B z < \K e ff\ because 



F = -M Z B Z cos 9 - K 2 cos 2 9 

where M z is the z-component of the magnetization and 
K 2 is the first nonvanishing coefficient in an expansion of 
the free energy for a system with second order anisotropy. 
For the equilibrium angle one gets: 



d 2 F{0) 
d6 2 



< for B z < \K eff \ 
> otherwise 



(A2) 



dF(9) 
88 



M Z B Z sind + 2K 2 cos 8 sin 8 = 0. (Al) 



holds. For a detailed discussion of stability conditions in 
film systems we refer to Refs. [i1l27l 
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